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We study S- wave resonance contributions to the —» J/ipTr + 'K~ and B s —» n + n~p + pT decays 

in the perturbative QCD framework by introducing two-hadron distribution amplitudes for final 
states. The Breit—Wigner formula for the /o(500), /o(1500), and /o(1790) resonances and the 
Flatte model for the /o(980) resonance are adopted to parametrize the timelike scalar form factors 
in the two-pion distribution amplitudes, which include both resonant and nonresonant contributions. 

The resultant branching fraction and differential branching fraction in the pion-pair invariant mass 
for each resonance channel are consistent with experimental data. The determined S'-wave two-pion 
distribution amplitudes, containing the information of both resonant and nonresonant rescattering 
phases, can be employed to predict direct CP asymmetries of other three-body hadronic B meson 
decays in various localized regions of two-pion phase space. 

PACS numbers: 13.20.He, 13.25.Hw, 13.30.Eg 


I. INTRODUCTION 

After the LHCb Collaboration measured sizable direct CP asymmetries in localized regions of phase space M- 
three-body hadronic B meson decays have attracted a lot of attention recently. To calculate direct CP asymmetries 
in the kinematic region with two final-state hadrons almost collimating to each other, we have proposed a theoretical 
framework based on the perturbative QCD (PQCD) approach with the crucial nonperturbative input of two-hadron 
distribution amplitudes [3, 6 ]. The typical PQCD factorization formula for aB-> /11/12/13 decay amplitude is written 
as @,|i| 


A = (j> B ® H ® (j) hlh2 ® <j > ha , ( 1 ) 

where the leading-order hard kernel H contains only one hard gluon as in the formalism for two-body B meson decays, 
and the B meson (/11-/12 pair, /13 meson) distribution amplitude 4 >b { 4 > h 1 h 2 i Sh 3 ) absorbs nonperturbative dynamics 
in the process. The analysis of three-body hadronic B meson decays is then reduced to that of two-body ones. 
Fitting the timelike form factors and rescattering phases contained in the two-pion distribution amplitudes to relevant 
experimental data, we have been able to make predictions for the direct CP asymmetries in the B± —> and 

7r + 7r _ Jv ± modes, which are consistent with the LHCb data in the localized region of the pion-pair invariant mass 
squared m 2 + w _ low < 0.4 GeV 2 Q. 

A limitation of the study in Ref. [ 7 | is that only nonresonant contributions to the timelike form factors were 
included, so it cannot be extended to regions involving intermediate resonances. In this paper we will propose 
parameterizations of the complex timelike form factors, which include both resonant and nonresonant contributions, 
taking the S'-wave two-pion distribution amplitudes as an example. The LHCb data of B s meson decays through 
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5-wave resonances, such as B s — > <//V , /o(980)[/o(980) — > 7r + 7r _ ] and B s — > f 0 (980)]/o(980) — > tt + tt ~}p + p~, are then 
adopted to fix involved parameters. It will be demonstrated that the resultant branching fraction and differential 
branching fraction in the piompair invariant mass for each 5-wave resonance channel agree with the data. It implies 
that the PQCD approach [8f 10] is an appropriate framework for analyzing three-body hadronic B meson decays. The 
determined 5-wave two-pion distribution amplitudes, containing the information of both resonant and nonresonant 
rescattering phases, can be employed to predict direct CP asymmetries of other three-body hadronic B meson decays 
like B± —> K ± 'k + 'k~ and B± —> n ± TT + n~ in various localized regions of two-pion phase space. The same formalism 
will be applied to the extraction of the P -wave two-pion distribution amplitudes from, for instance, the Belle data 0 
in the future. 

The B s —> J/ipTr + TT~ decay was first measured by the LHCb Collaboration [l2j], one of whose channels B s —> 
J/^/>/o(980)[/o(980) —> 7r + 7r - ] with a CP-odd final eigenstate can be used to extract the mixing-induced CP violation 
phase (j) s . A modified Dalitz-plot analysis [l3| indicated that the /o(980) and /o(1370) resonances give the dominant 
intermediate contributions, about 70% and 21% of the total decay fraction, respectively. It was shown based on 
3 fb -1 of integrated luminosity 0 that the 5-wave resonances /o(980), /o(1500) and a so-called /o(1790) state, 
along with tiny components from the D -wave resonances /2(1270) and /2(1525), well describe the B s —> J/ipn + TT~ 
data. The Cabibbo- and color-suppressed decay B° —> J/ipn + TT~ was first observed by the BABAR Collaboration 
with the branching ratio (4.6 ± 0.7 ± 0.6) x 10” 5 , to which the B° —> J/ipp 0 channel was found to contribute 
(1.6±0.6±0.4) x 10 -5 [Ij|, and (2.7±0.3±0.17) x 10 -5 in a later improved measurement 0]. The LHCb Collaboration 
confirmed these results by investigating resonant components in the B° —> J/ipTr + n~ mode thoroughly 00: the 
intermediate p° meson dominates, contributing about 63% of the total decay fraction EH , and the 5-wave resonance 
/o(500) is the next, giving about 22% EH- In contrast to the sizable /o(500) component, another 5-wave state /o(980) 
was not seen in the B° J/tpTr + TT~ decay 00- 

The first observation of the B s —> tt + tt~ p + p~ decay and the first evidence of the B° —> p + p,~ decay have 

been announced by the LHCb Collaboration 0], with the branching ratios (8.6 ± 1.5 ± 0.7 ± 0.7) x 10“ 8 and 
(2.11 ± 0.51 ± 0.15 ± 0.16) x 1CD 8 , respectively. Within the measured range of the pion-pair invariant mass 0.5- 
1.3 GeV, the B s —> / 0 ( 980)p + p~ ( B° —► p(770)°p + p~) channel is expected to dominate in the former (latter), 
contributing (8.3 ± 1.7) x 10 -8 ((1.98 ±0.53) x 10 -8 ). Hence, the B s —► TT + n~p + p~ data also provide a useful input 
for the determination of the 5-wave two-pion distribution amplitudes. It will be shown that they are crucial for fixing 
the involved Gegenbauer coefficient, to which the B s —> data are less sensitive. For an independent study of 

the B° ir + ir~p, + p,~ decay in the light-cone sum rules, refer to Ref. EH- The B° p(770) o [p(770) 0 —> TT + n~]p + p,~ 
decay is applicable to the extraction of the P- wave two-pion distribution amplitudes stated before. 

The factorization formalism of the B®^ —> J/'i/>7r + 7r _ and B s —> p + p,~ decays with the inputs of the 5-wave 

two-pion distribution amplitudes is established in Sec. II. The associated timelike scalar form factors are parameterized 
in terms of the Breit-Wigner (BW) model for the /o(500), /o(1500) and /o(1790) resonances and the Flatte model for 
the /o(980) resonance. In Sec. Ill the relative strengths and strong phases among the above resonances are determined 
by fitting the factorization formulas to the LHCb data. It will be highlighted that the obtained differential branching 
fraction in the pion-pair invariant mass for each resonance channel matches well the experimental one. Especially, our 
results for the B s —> J/-i/>7r + 7r _ decay favor the ’’Solution I” set of branching fractions for various resonances presented 
by the LHCb Collaboration 0- Section IV contains the conclusion. The explicit PQCD factorization formulas for 
all the decay amplitudes are collected in the Appendix. 


II. DECAY AMPLITUDES 

A. B%) 7/V’ 7r+7r_ Decay 

The B® s j —> J/il)TT + n~ modes mainly proceed via quasi-two-body channels containing scalar or vector resonant states 
as argued in Refs. S0- For discussions on the B s —> J/tj) 

The chiral unitary approach including final state interaction 
decays recently. The B° —> branching ratio was calculated in the QCD-improved factorization approach, 

where a two-meson distribution amplitude for the pion pair and final state interaction were considered [29| . In this 
section we derive the PQCD factorization formulas for the B^ —t decays with the inputs of the 5-wave 

two-pion distribution amplitudes. It has been postulated that the leading-order hard kernel for three-body B meson 
decays contains only one hard gluon exchange as depicted in Fig.[l] where the B° or B s meson transits into a pair of the 
7 r + and 7r _ mesons through an intermediate resonance. Figures []Ja) andQJb) represent the factorizable contribution, 
and Figs. [Qc) andQJd) represent the spectator contribution. 

In the light-cone coordinates, the B ? s ^ meson momentum ps, the total momentum of the two pions, p = p\ + P 2 , 


channel, refer to Refs. [21-125 
was employed to investigate the B^ J/ip7r + TT 
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JW JN> 



(a) (b) 


JM> J/V 




FIG. 1: Typical Feynman diagrams for the three-body decays B®^ —> J/ipn+n , where B stands for the B° or B a meson, and 
® denotes the weak vertex. 


and the J/ip momentum p 3 are chosen as 

m B 


VB = ~^ 1,1,0t )’ P = “ r ’ 2 ’ 7 ?’° T )’ P 3 = ^|(r 2 ,l - 77,0 T ), 

where m B denotes the 7 >'F meson mass, mj is the J /ip meson mass, and the variable 77 is defined as 


( 2 ) 


OJ 


77 = 


(1 — r 2 )m 2 B 


2 > 


(3) 


with the mass ratio r = mj/^/m B and the invariant mass squared w 2 = p 2 of the pion pair. Define C, = p+ /p + as the 
7 r + meson momentum fraction, in terms of which the other kinematic variables of the two pions are expressed as 


Pi = (! - pi = P2=c r i'—§t- 


V2' 


m B 


(4) 


The momenta for the spectators in the B ^ meson, the pion pair, and the J /ip meson read as 

'm B 


k b = (0 ,—j=x B ,k B T 


k = ( ~ r2 )’°> fc T 


k 3 = (vf r2x3 ’ ^ _ r, ' )X3,k3T 

where the momentum fractions x B , z, and x 3 run between zero and unity. 
We introduce the distribution amplitudes for the pion pair 


(5) 


$w(2hC,^ 2 ) = 

*te,C,w 2 ) = 


1 


d V 


2 v / 2 W c 

J- 

2tt e 

1 

P+ 

f dy ~ e- 

2y/2W c 

w 

J 27T 

1 

ft 

[ d y~ e 

2 V 2 w c 

w 2 

j 2 7T 


'Hi/' 


( 6 ) 

(7) 

( 8 ) 


where T = 1/2 (T = r 3 /2) is for the isoscalar 7 = 0 (isovector 7=1) state, N c the number of colors, ip the u-d doublet, 
z the momentum fraction carried by the spectator u quark, and n_ = (0,1, Oy) a dimensionless vector. For 7=1, 


11 ir: j. -wcivc 10 11 Lr- irau 111 |jai iiai wave, anu 'm_ auu econtribute at twist-2. The other components 4*(, 
<Ff =1 , and then contribute at twist-3. The v = +, —, and T components can be extracted by contracting the 

above definitions with p + , n_, and pi Bl respectively. For 7 = 0, the S-wave is the leading partial wave, and $ 7 / 7 £_ 
contributes at twist-2. The others and contribute at twist-3. In the case with the pion pair coming from 

the s and s quarks, ip is replaced by the s quark field. A two-pion distribution amplitude can be related to the pion 
distribution amplitude through the evaluation of the process 77 * —> at large invariant mass w 2 [Hil. The 

extraction of the two-pion distribution amplitudes from the B —» ttttW decay has been elaborated in Refs. pT PG|. 

The 5-wave two-pion distribution amplitudes are collected into 37] 


$ 


S—wave 




i^L(z,C,w 2 ) + w $f=°( 2 ,c, w 2 ) + - 1)C=° + (*,C,™ 2 )], 


71+ = (1,0, Or) being a dimensionless vector. Their asymptotic forms are parameterized as 


$ 




${=° = 


F s (w 2 ) 

2y/2N~c 


T, /=0 , 


F s (w 2 ) 

2^2W c 


( 1-2 z), 


(9) 


( 10 ) 























4 


with the timelike scalar form factor F s (w 2 ) and the Gegenbauer coefficient a^ _0 . Following the analysis of the LHCb 
Collaboration 0 BEE 0, and motivated by the studies in Refs. [38l - i40j . we introduce the S-wave resonances 
into the parametrization of F s (w 2 ), so that both resonant and nonresonant contributions are included into the S'- 
wave two-pion distribution amplitudes. For the ss component in the B s —► J/ipT: + TT~ decay, we take into account 
/o(980), /o(1500), and /o(1790) as in Ref . 0 . It was noticed by the LHCb Collaboration that the above resonant 
states lead to a fit better than in Ref. 0 , where /o(980), /o(1370), and /o(1500) were considered. Note that there is 
a resonant state /o(1710), instead of /o(1790), in Refs. [4l|, |42(] with the mass around 1.72 GeV. 

For the / o (1500) and /o(1790) contributions, the BW formula is employed 


m| — to 2 — imsTs(uj 2 ) ’ 

where ms is the pole mass of the resonance, and the energy-dependent width r,g(u; 2 ) for a S-wave resonance decaying 
into two pions is parameterized as 


r s (w 2 ) = r s 


w 2 — 4m 2 \ 
w V m s — 4m 2 ) 


F-, 


Ri 


( 12 ) 


with the pion mass m n , the constant width Tg, and the Blatt-Weisskopf barrier factor Fr = 1 in this case [13j, |43|. 
The values of m/ 0 ( 1500 ) = 1-50 GeV and T / 0 (i 5 oo) = 0-12 GeV are taken for the mass and the width of the /o(1500) 
meson, and toj 0 ( 1790 ) = 1.81 GeV and T j o ( 1790 ) = 0.32 GeV are taken for / 0 (1790) 0, El], |42|] . The BW formula does 
not work well for / 0 (980), since this meson is close to the I\ K threshold. We then employ the Flatte model [44j , 


m fo(980) ^ ^^/oCSSO) {SnirPTTK T 9KKPKk) 


(13) 


with the pole mass my o ( 980 ) = 0.97 GeV [41|, |42|, |45|, |46|. The coupling constants g nT7r = 0.167 GeV and gKK = 
3.47 5 ™ m 0 describe the /n(980) decay into the final states 7r + 7r and K + K , respectively. The phase space 
factors p 1Tn and pkk read as [lj, 0,01 


Ptx 


2 / 4m 2 ± 

= -\ 1 -^ 


U 




cu* 


1 L 4 m 2 K± 
PKK = 2\ 1 -^^ 


y 1 - im 


UJ 


(14) 


Considering the relative strengths and strong phases among different resonances, we write the timelike scalar form 
factor associated with the ss component as 


*7K) = 


cim 2 o(980) e iei 


W /o(980) w2 * TO /o(980) {d^TTP-ir-K + gKKPKK) 


+ 


^^(ISOO )^ 2 


m /o(1500) ~ w2 - *"l/o(1500)r /o (l500)(w 2 ) 


+ - 


C 3 m /o(1790) e 


m 


/o(1790) 


- w 2 - nn /o( 1790 )r /o( 1790 ) (u; 2 ) ’ 


(15) 


Ci and 9i, i = 1, 2, and 3, being tunable parameters. The above parametrization decreases like ui~ 2 in the asymptotic 
limit, indicating that it also includes the nonresonant contribution Q. For the dd component in the S'-wave two-pion 
distribution amplitude, only the resonance /o(500) or the so-called a meson in literature is needed. It was observed 
that the /o(980) contribution to the_R° —> J/^/> 7 r + 7 r - branching ratio is less than 1% compared with the /o(500) 
one E3- That is, the LHCb data HzL0 disfavor the existence of /o(980) in the B° —> J/i/jtt+tt decay, rendering 
doubtful the hypothesis of a four-quark state (with the dd content) for this meson. The /o(500) contribution to the 
scalar form factor for the dd component is parameterized as a BW model with the pole mass m/ 0 ( 500 ) = 0.50 GeV 
and the width T / 0 ( 5 oo) = 0.40 GeV according to Refs. [13 [H, Hi], 0,0, 


F. fV) 


_ CW /o(500)_ 

m /o(500) ~ 'F 1 ~ * m /o(500)r/o(500)( w2 ) 


(16) 


c being a tunable parameter. 
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The differential branching ratio for the B H — > J/ipir + n decay is expressed as fill. [~12l 

^ _ ^IpiIIpsI i a 1 2 
duj B 4(2n) 3 m B ’ 


(17) 


with the meson mean lifetime t b ■ The kinematic variables \pt\ and |p|| denote the magnitudes of the 7 r+ and 
J/i/> momenta in the center-of-mass frame of the pion pair, 



4m 2 + , 

•7T n= 7 


1^1 



(u + mj/^) 2 ] [m| 


(w - m J/v ,) 2 ]. 


(18) 


The decay amplitudes A from Fig. [T]with the 5-wave pion pair in the final state are collected in the Appendix. 


B. B s -A 7r + 7r p+p Decay 


Next we establish the PQCD factorization formalism for the differential branching ratio of the B s —> tt+7t p, + p 
decay, which is generically a four-body mode. This decay is governed by the effective Hamiltonian 


Heft = ~^v tb vz CiMOiM, 


10 


(19) 


where Oi and C, are a four-quark or magnetic-penguin operator and its corresponding Wilson coefficient, respectively, 
the renormalization scale p will be set to the b quark mass mb below, Gf is the Fermi constant, and Vtb and Vt s are 
the CKM matrix elements. Equation (fTTFl) yields the effective operators for the b -A s£ + £~ transition 


iM{b ~A s£+r) = iNA (Cg + C 10 )[sb\v-A[U]v+A + ( Cg - C W )[sb\ V . A [U] V - A 




-\-4C 7L m b -^[sia^(l + TsM^y^] + 4C 7R m b -^[sia ilv (l - 7 5 ) 6 ][Hy I '£] 
with C 7 l = C 7 , C 7R = C 7 Lm s /mb , m s being the s quark mass, the lepton-pair momentum q , and the coefficient 

m = %^v tb v t * s . 


( 20 ) 


4\/2 


( 21 ) 


The B s -A p + p~ decay amplitude is obtained by sandwiching Eq. (1201) between the initial and final hadronic 
states, to which only the factorizable diagrams, Fi gs, pl a) andQJb), contribute. The spinor products of s and b are 
then replaced by the hadronic matrix elements [481149j | 


((tt+tt )s(p)|s7m75&|H s (pbJ) = ~i 


P»- 


TOr — U) 


-9m 


((tt+tt )s{p)\s(T^q v l5b\B s {pB a )) = 


(m| - u 2 )q fJ , - q 2 P^ 


Pi(q 2 ) + 

Hr (q 2 ) 

. u(iti Bb + w) ’ 


Trim — UJ 


-q^oiq 2 ) c 


( 22 ) 


that define the B s -A tt+7t~ form factors T o, Hi, and J~t with the sum of the momenta P = p Ba + p. The 
explicit factorization formulas for the above form factors are referred to the Appendix. Since the PQCD approach is 
applicable only to the large recoil region, the extrapolation of the form factor behavior to small recoil is necessary for 
the evaluation of the branching ratio. Here, we adopt the dipole form for the the q 2 dependence of the form factor 
Hi, i = 0, 1, and T, 


Pi(q 2 ) = 


Hi(0) 


1 - a l g 2 /m| + bi(q 2 /m 2 B ) 2 ’ 


(23) 


where the shape parameters cq and 6j, depending on a; 2 , will be obtained from the knowledge in the range q 2 < lOGeV 2 . 

A general differential decay width for the B -A h\h- 2 (- + b~ mode with various partial wave contributions has been 
derived using the helicity amplitudes in Refs. [50l.[5l|. In the case with the dominant 5-wave contribution, the angular 
distribution is written as 


d 3 r 


3 r 


duj 2 dq 2 d cos 0e 8 L 


+ J 2 C cos( 2 ^) 


(24) 
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with Oe being the polar angle between the £ and B s moving directions in the lepton-pair rest frame. The angular 
coefficients are given by 


Jf = 


l^o! 2 + l^ol 


SmA-AloA'gol cos {6% - 5° m ) + 4m e \A" t 


0 I 2 


JH = -t 


1-4' 


L Ol 


1-4 


R 0 \ 


(25) 


with Pi = \J\ — 4m'j/q 2 and me = me/^/q 2 , in which the subscript t denotes the timelike component of a virtual 
state decaying into a lepton pair, and the strong phases S B0 and <5^ 0 associated with the helicity amplitudes Alo and 
Arq, respectively, vanish at leading order in the strong coupling a s . The helicity amplitudes are expressed, in terms 
of the form factors, as 


A 0 

L/R ,0 


= V^2- 


Vx 


V\m Bs 


(C 9 q= + 2(C 7L - Ctr) _ 

q 2 VQ ( m Bs+ U} ) 


A T (q 2 


A° t = 2 y/W 2 C w - 

UJ 


V? 


-Mo 2 


(26) 


with the factor 


No = 


167T 2 


NiNnnVl - 4m 2 /w 2 , = 


8 y/\q 2 


3 2567r 3 m^ a 


(27) 


and the function A = A(m| , w,q 2 ) = {m 2 B — w 2 — q 2 ) 2 — 4 uj 2 q 2 . 


III. NUMERICAL RESULTS 


In addition to the quantities that have been specified before, we take the following inputs (in units of GeV) i7. FTl 

A^Z 4) = 0.250, m B0 = 5.280, m Ba = 5.367, m J/v , = 3.097, 
m,± = 0.140, 771^0=0.135, nib = 4.66, m s = 0.095, (28) 

the mean lifetimes r R o = 1.519 x 10 -12 s and t Bs = 1.512 x 10 -12 s [U, |42[, and the values of the Wolfenstein 
parameters in Refs. |4jJ,|42|. Choosing the Gegenbauer coefficient as = 0.2, we extract the parameters from the 
LHCb data 0, 


ci = 0.900, c 2 = 0.106, c 3 = 0.066, c = 3.500, 

0i = 0 2 = 03 = 0 , 


which correspond to the branching ratios 

B(B S J/t/>/ o (980)[/o(980) ->■ tt+tt-]) 
B{B S -a J/V’/o(1500)[/o(1500) -»■ 7 t+tT]) 
B(B S -»■ J/^/ o (1790)[/o(1790) -)■ tt+tt-]) 
B(S° J/V>/o(500)[/ o (500) 7r+Tr"]) 


(i-is±8;S 

(l-62±g-iS 

(3.26™ 

(6.91±1;S 


(«B.) 18 :ig(^=°)l 8 :S(mc)) x io- 4 , 

(^Jt 0 o-J 8 o(«r°)iro^(mc)) X 10- 5 , 
(^.n.SS(^=°)l8:8S(mc)) x io- 6 , 
(«B)liS(^= 0 ) 18 :g(»nc)) x io- 6 , 


(29) 


(30) 


with the three errors coming from the variations of uis a = (0.50 ± 0.05) GeV [52j or ui B = (0.40 ± 0.04) GeV, 
a 2 ~° = 0.2 ± 0.2, and m c = (1.275 ± 0.025) GeV. It is seen that the above results are not sensitive to a 2 =0 , so 
the central value of a 2 ~° is actually mainly determined by the B s —» data. The errors induced by the 

variations of the Wolfenstein parameters and the mean lifetime of the B° or B s meson are tiny, and have been omitted. 

We may set the B s —> J/t/j/o(1500)[/o(1500) -A 7r + 7r~] contribution to be 10.0% of the total B s —> J/i/m + TT~ decay 
rate, which roughly agrees with the corresponding values in both Solutions I and II presented by the LHCb Collabo¬ 
ration 0. The percentages of the other channels are then given by 70.9% and 2.0% for B s —> J/t/j/o(980)[/o(980) — > 
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7r + 7r ] and B s — > <//'i/’/o(1790)[/o(1790) — > tt + tt ], respectively, and support the Solution I data [HJ- Including all 
the 5-wave resonances / o (980), / o (1500) and / 0 (1790) in the scalar form factor, we have the total branching ratio 

B(B S -+ J/^(n+n-) s ) = (l.62l“;“(u;B.)l“;“( O r 0 )ig:^(m c )) x 10" 4 , (31) 

to which the interference between /o(980) and /o(1500) (/o(980) and /o(1790), /o(1500) and /o(1790)) contributes 
1.61 x 10 -5 (1.42 x 10 -5 , —2.13 x 10 -6 ). The B° —>■ J/^/ o (500)[/ o (500) — > 7r + 7r~J branching ratio predicted in Eq. (1501) 
is consistent with the data (8.8 ± 0.5+g) x 10~ 6 in Ref. [17) and (6.4 ± 0.8+ s ) x 10 -6 in Ref. )18j |. 

For comparisons with other data, we read the branching ratio B(B S —> J/ipn + TT~) out of the ratio [l3| 

= (19.79 ±0.47 ±0.52)%, (32, 

with the branching ratio B(B S —> J/ip(j>) = (1.07±0.09) x 10~ 3 [ill.l42j|. Then a value B{B S -+ J/i/jf o (980)[fo(980) —> 
7r+7r_ ]) = (1.48 ± 0.14) x 10 -4 is inferred from its percentage in the total decay rate of B s J/i/>7r + 7r _ [13[. This 
outcome and the Belle data B(B S -+ </+/o(980)[/o(980) — > 7r + 7r - ]) = (1.16+ 3 g(stof)lo+(s2/sf)lo as(N B (.») )) x 
10 -4 [H[ are both compatible with the PQCD prediction in Eq. OTTj) . 

A ratio of the decay rates, 


□ = r (B s ° 


J/lpfoJo 7T+7T ) 
T(B° -> <j> -+ K+K-) 


= 0.252 


+0.046+0.027 

-0.032-0.033) 


with /o standing for the / 0 (980) meson, was measured by the LHCb Collaboration 
by the CDF and DO Collaborations [H, [13 ■ We predict this ratio 


td _ n 99 /A+ 0 .094+0.034+0.004+0.020 

n fo/4> u - z - jU -0.073-0.029-0.002-0.017) 


(33) 

|, and subsequently confirmed 

(34) 


where the first (second, third) error comes from the variation of the shape parameter ujb s (the Gegenbauer moment 



(a) (b) 

FIG. 2: Dependencies of the differential branching ratios dB/dui on the pion-pair invariant mass for (a) the resonances /o(980), 
/o(1500), and /o(1790) in the B s -+ J/ip 7r + 7r _ decay, and (b) /o(500) in the B° -+ J/?/,7r + 7r _ decay. 


The predicted dependencies of the differential branching ratios dB/duj on the pion-pair invariant mass oj are pre¬ 
sented in Fig. [2ja) for the resonances /o(980), /o(1500), and /o(1790) in the B s -+ J/ilm + Ti~ decay, and in Fig. OJb) 
for /o(500) in the B° -+ J/ipn + n~ decay, which can be compared with the data in Refs. [Til . Til . IT . Il8 : ] . A value of 
the differential branching ratio at each invariant mass of the pion pair is proportional to the event number observed 
by the LHCb Collaboration. We have adjusted the range of the y axis, such that it is easy to see the coincidence of 
the the curves in Fig. [2] and the corresponding ones in Refs. (ill. ITtI . 

The integrated branching fraction for the B s — > n + tt~ fi + y,~ decay was measured by the LHCb as [fj| 

B{B S -+ 7T+7 r-y+fj,-) = (8.6 ± 1.5 ± 0.7 ± 0.7) x 10" 8 , 


(35) 














FIG. 3: (a) Predicted differential branching ratio for the B s —» 7r + 7r ~[x + decay, (b) Comparison between the experimental 
data (with triangle markers) the theoretical predictions (with square markers) for the differential branching ratio of the B s —► 
7r + 7r - n~ decay. 


where the first two errors are statistical and systematic, respectively, and the third error is due to the uncertainty 
from the normalization, i.e. the B° —> -± p, + /i _ )A '*(—> K + ir~) branching fraction. Adopting the parameters of 

the 5-wave two-pion distribution amplitudes fixed in Eq. (l29l) . we calculate Figs. [Ha) andQJb), and the differential 
branching fraction for the B s —> 7r + n~ p, + /i _ decay. Our result is displayed in Fig. E[a), which clearly exhibits the 
peak arising from the /o(980) resonance. Integrating over ui, we derive 

B{B a TT+n-fi+n-) = (8.0 ± 2.4(w s J ± 2.5(4=°)) x 10“ 8 , (36) 

consistent with the data in Ea. (l35l) . where the theoretical errors are attributed to the variations of u>b b = (0.50±0.05) 
GeV and 4 = ° = 0.2 ± 0.2. In order to be compared with the LHCb data [lj$] carefully, we also present the binned 
results in Fig. [3])b) from 0.5 GeV to 1.3 GeV, in which the data (with triangle markers) have been normalized to the 
central value B(B° S —► 7r + 7r - /i + /i - ) = 8.6 x 10 -8 . The comparison shows a general agreement between our theoretical 
predictions and the data. The different around w ~ 0.8 GeV, if becoming more significant in the future, may imply the 
inclusion of the missing /o(600) resonance into the parameterized form factor. It is noticed that the B s —> 7r + 7r - /i + /i - 
mode is more sensitive to the parameter 4=° than B s —> J/i/jTT + Tr ~, so the comparison in Fig. [^b) is more crucial 
for the determination of 4=° as stated before. 

With the assumption on the dominance of the /o(980) —> transition, the B s —> / o (980)/i + /r“ branching 

fraction was found to be 0 

B(B S / 0 (980)(—i tt + tt-)h + ii-) = (8.3 ± 1.7) x 1CT 8 . (37) 

On the other hand, the branching fraction B(B S —> J/ipf 0 ) = (1.39 ± 0.14) x 10~ 4 (42| would indicate 

B(B S /o(980)/r"4 _ ) ~ 6.1 x lO" 8 , (38) 

according to the ratio [42} 


B{B~ K~n+ir) (4.49 ± 0.23) x 10" 7 
B{B~ -)• J/4’K~) ~ (1.027 ±0.031) x 10“ 3 


(39) 


if the above ratio was not sensitive to light final-state mesons as in most cases. The result in Eq. (1381) is consistent 
with that in Eq. J37l), but the central value is lower by about 30%. Both Eqs. (l37l) and (l38l) agree with our prediction 
in Eq. (1M1) basically. More precise data from the LHCb and the future KEKB will test our theoretical formalism. 
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IV. CONCLUSION 

In this work we have studied the contributions from the 5-wave resonant states /o(980), /o(1500) and /o(1790) 
to the B s —> decay, and from / 0 (500) to the B° —> J/iI’tt + tt~ decay in the PQCD approach with the 

introduction of the S'-wave two-pion distribution amplitudes. The BW formula for / o (500), / o (1500) and / o (1790) 
and the Flatte model for /o(980) were adopted to describe the involved timelike scalar form factors, which contain 
both resonant and nonresonant dynamics in the two-pion system. The resultant branching fractions, ratio Rfg/^, and 
differential branching fractions in the pion-pair invariant mass for the relevant resonance channels all agree with the 
data well. The same two-pion distribution amplitude associated with the ss component also leads to the differential 
branching fraction of the B s —> decay, which is compatible with the data. It is stressed that the data 

of both B s —> J/7t~ and B s —> 7r + n~ /i _ modes are required for determining the timelike scalar form factor 
and the Gegenbauer coefficient in the S -wave two-pion distribution amplitudes completely. We conclude that our 
formalism with the two-hadron distribution amplitudes is successful, and applicable to other three-body hadronic 
B meson decays. The extracted S -wave two-pion distribution amplitudes can be employed to predict direct CP 
asymmetries of the decays with intermediate S -wave resonances in various localized regions of the two-pion phase 
space. 
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Appendix A: Decay amplitudes 


In this Appendix we present the PQCD factorization formulas for the diagrams in Fig. [T) The J/ij) vector meson 
distribution amplitudes with longitudinal polarization read as 


* L = 1 

3 ^ 




(Al) 


with the longitudinal polarization vector 


£l = 


m B 


\/2 rrij/^ 


{—r 2 , (1 — rj), 0 T ) • 


(A2) 


The twist-2 distribution amplitude ip L and the twist-3 distribution amplitude ip* take the forms 


ip L (x) = 9.58 — x ) 

r V ' 2y/2Nc y ’ 

= 10.94 f J /^L {l-2x) 2 
y K 1 2 v / 2W 


x(l — x) 


1 — 2.8x(l — x)_ 


c(l — x) 


0.7 


n 0.7 


1 - 2.8x(l - x) 


(A3) 


with the J 1 1 \) meson decay constant fj/^. The distribution amplitudes and the relevant parameters for the B° and 
B s mesons can be found in Ref. ,58|. 

The analytic formula for the B —> J/^/)7r + 7r _ decay amplitude is given by 

A(B° [S) = V; b V cd(cs) ( F ll + M ll ) - V t * b V td{ts) ( F' LL + F LR + M' LL + M SP ) . (A4) 

For simplicity, we denote the distribution amplitudes (z, (, w 2 ) [*3>f ==0 ( 2 ;, Q w 2 ), w 2 )] by </>o (4>s,<t>a) 

below. The amplitudes for the B ( B° or B s ) meson transition into two pions from Fig. [Ija) and Fig.[ljb) are written 
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as 


p 1 poo 

F LL = 8nC F m 4 B fji^ / dx B dz / b B db B bdb(j) B (x B ,b B ) 

JO Jo 

x | VvO - r 2 ) [((1 - 2z)(l -rj)- r 2 ( 1 - 2z(l - ??)))(<^ s + <j> a ) + 2r 2 </> CT ] 


+ [(1 + z)( 1 - 77) - r 2 (l + 2z(l - ??))] 

+ 


ai(ta)E e (t a )h a {xB, z, b B , &) 


2 \A?(! - r 2 ) [1 - r) - r 2 ( 1 - x B )] 
-(! - r 2 ) [(1 - ??)f 7 + r 2 {x B ~ v)] 0 o 

p/LL _ p' ' 
j^LR _ ^LL 


ai(t b )E e (tb)h b (x B , z, b B , b ) 


I Cti —>-0.2 > 


(A 5 ) 

(A6) 

(A 7 ) 

with the Wilson coefficients ai = C\ + C2/N c , a,2 = C 3 + Ci/N c + Cg + C\o/N c , and 03 = C5 + Cq/N c + C7 + C$/N c . 
The spectator diagrams in Fig. OJc) and Fig. [ljd) yield the amplitudes 


a\—^a^ ■> 


m ll = - 


pi poo 

32fCfbi%/VG / dx B dzdx 3 / b B db B b 3 db 3 (f> B (x B ,b B ) 

Jo Jo 

VvO - r 2 )[2(l - 77) (z( 1 - r 2 ) + r 2 (l - £3))^ 

-(z(l - 77) (1 - r 2 ) + r 2 x s )(0 s + 0 < t)]' 0 L + (1 - r 2 - 77) 
x [ ((1 £3 Xb)( 1 - r 2 ) + 77(2:3(1 - 2r 2 ) - (1 - r 2 )(l - z) + r 2 )) ip L 


+rr c { 1 - 77)^] 0o 


C2{t c )E n {t c )h c (x B ,z 1 x 3 ,b B ,b 3 ) 


- 77(1 - r 2 )[( 2 r 2 (xs - £3(1 - + 4 rr c (l - r])^)^ 

-{r 2 x B + z( 1 - 77)(1 - r 2 )){(j) s + 0 (T )0 i ] 

-(1 - r 2 - 77) [(x B - z - x 3 (l - 77) + r 2 (z - x 3 ))tp L + rr c ( 1 - 77)7/0] 0o 

xC2{td)E n (t d )h d {x B ,z, x 3 , b B , b 3 ) 

M' LL = M LL \c 2 ^ ai , 

pi poo 

M SP = 32 nCFm B /V 6 / dx B dzdx 3 / b B db B b 3 db 3 (j) B (x B ,b B ) 

Jo Jo 

- \J 77(1 - r 2 )[ 2 r 2 ((l - x 3 )(l - 77) - x B )il> L - 4 rr c (l - ?7)0*]0 O 


+ ^77(1 - r 2 ) [z(r 2 - 1)(1 - 77) - r 2 £ S ] (0 S + 0<r)7/’ i: ' 

+(1 - r 2 - 77) [((1 - x 3 )(l + r 2 - 77) + z(l - r 2 ) - x b )i/j l 


—rr c (l - 77 )' 0 t ] 0 o 


^5 (tc)E n (Jc)hc(^x B , z, £3, b B , 63) 


>/7?(1 - r 2 ) [2(77 - l)(r 2 (x 3 - z) + z)<j> a + (jz(1 - r 2 )(l - 77) + r 2 x B )(<fr s + <j>„)\ i/j l 
+(1 - r 2 - 77) [ ((1 - r 2 ){x B - zrj) + x 3 (r 2 (l - 2rj) - (1 - 77))) i/j l 


—rr c (l - 77)^] 0o 


a 5 (td)E n (td)hd{x B ,z, x 3 , b B ,b 3 ) 


(A8) 

(A 9 ) 


(A 10 ) 


with the mass ratio r c = m c /m B , the notations = ^ L,t (x 3 ), and the Wilson coefficients = C4 + C w and 
®5 = Cq + C 8 . 
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The hard functions h a ,hb, h c and h d are written as 


h a (x B ,z,b B ,b) 


hb(x B ,z, b B , b ) 


h c (x B ,z,x 3 ,b B ,b 3 ) 


h d (x B ,z,x 3 ,b B ,b 3 ) 


[9(b- b B )I 0 (m B s/z(l - r 2 )b B )K 0 (m B ^/z{l - r 2 )b) + (6 6 B )] 

xKo(m B y/ab B )S t (z), 

f [0(b-b B )I 0 (m B ^/Kb B )K 0 (m B ^/Kb) + (b <-> b B )\, k>0 


ITT 

2 


9(b - 6 B ) J 0 (m B v^6 b )F^ 1) (to b y / j/^6) + (b b B ) 


k < 0 


x K 3 {vn B y/ob^jSt{\x B - rj\), 

[0(b B - b 3 )I 0 {m B y/ab 3 )K 0 (m B y/ab B ) + (6 B -n- 6 3 )] 
( K 0 (m B y/]5b 3 ), (3> 0 

(3 < 0 


h c (x B ,z, 1 - x 3 ,6 B ,6 3 ), 


(All) 


with the factors a = (1 — r 2 )x B z , k = (1 — r 2 )(x B — ij), and /3 = r 2 — (z( 1 — ?' 2 ) + r 2 (l — x 3 ))((l — 77 )(1 — x 3 ) — x B ), 
and the Hankel function Hq 1 \x) = Jq(x ) + iYo(x). In the above expressions the threshold resummation factor St(x) 
follows the parametrization in Ref. [59( 


with the fitted parameter 


S t {x) 


2 1+2 T(3/2 + c) 

V / H 1 (l + c) 


[cc(l 


c = 0.04Q 2 — 0.51Q + 1.87, 

and Q 2 = m B ( 1 — r 2 ) or Q 2 = m 2 B { 1 — q 2 /m 2 B ). 

The evolution factors E e [t) and E n {t) in the factorization formulas are given by 


(A12) 


(A13) 


E e [t) = a s (t)exp[-S B (t) - S M (t)\, 

E n (t ) = a s (t)exp[-S B (t) - S M (t) - Sj/^(t)]\ b= b B , (A14) 

in which the Sudakov exponents are defined as 

Sb = s^XB^bs'j+lJ^ ^ 7 q (a s {fi)), (A15) 

s M = s ( z (i - f2 )^f ’ & ) + s ((i - 2 )(i - r2 )^f ’ 6 ) + 2 ( A16 ) 

Sj/4 = s (z 3 (l - ?y)^|-,6 3 ^ + s^(1-x 3 )(1-?7)^|,6 3 ^ + 2^ -|-7 ' q {a 8 {p)), (A17) 


with the quark anomalous dimension = — a s /Tr. The explicit expression of the function s can be found, for example, 
in Appendix A of Ref. tjfl]]. The hard scales t a ,tb,t c and t d involved in Figs. [Ha)-[Ud) are chosen as 


t c 


max 

max 


jm B \A(l - r 2 ), 1/6, l/6 B j , 

{ m B y/a , m B ^I,1/6 b ,1/6 3 }, 


tb = max jm B \/[ft[, 1/6,1/6 B | 


td — t c \x3—t(l—X3)- 


(A18) 


The PQCD factorization formulas for the B s —> tt + tt form factors with the S-wave two-pion distribution amplitudes 
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are written as 


/* 1 POO 

Fo(q 2 ) = 87 tCf^i 2 b s / dxsdz / bBdbBbdb(/)B s (xB^B) 

S Jo Jo 

x | [77(277 + 1)00 + 001 - r 2 ) r](l - 200 + 77(1 - r 2 )(- 2 zr ] - 77 + 2)0 S ] 

xh e (x B ,ZT], b B ,b)a s (tl)exp[-S B (tl) - ^JV^(*e)]‘ 5 'i(- 2r ) 

+ 20??(1 - r 2 )r]<t>sh e (z, x B rj, b, b B )ct s (t 2 e )exp[-S B (t 2 e ) - Sm^S/zb)} , 

/»1 /»oo 

J 7 !^ 2 ) = 87 tCftti 2 b s / dxsdz / bsdbsbdb(f>B s (^s, 6 b) 

s io ./o 

x | [07? + 1)00 - 1 / 77(1 - r 2 )(l + 2z- 2 / 77)0 - 077(1 - r 2 )(2z - 1)0] 
xh e (x B ,zr],b B ,b)a s (tl)exp[-S B (tl) - S M {tl)]S t (z) 

+20 77(1 - r 2 )(j) s h e {z, x B r], b, b B )a s {t 2 e )exp[-S B (t 2 ) - S M (tl)]S t (x B )} , 

Tr{q 2 ) = 8^607770 (1 + 077(1 - r 2 )) f dx B dz 

Jo 

x | [0?7(1 - r 2 )(-00 + 0 + 077(1 - r 2 )0 + 2 / 77 ) 0 ] 
xh e (x B ,zri 1 b B ,b)a s (t\)exp[-S B (t],) - S M (tl)]S t (z) 

+2077(1 - r 2 )cj) s h e (z, x B r], b, b B )a s (t 2 e )exp[-S B (t 2 e ) - Sm^)]^ (x b )| . 


b B db B bdb(j) Bs (x B ,b B ) 


(A19) 


(A20) 


(A21) 


The above expressions are similar to those for B meson decays into a scalar meson [62j except the replacement of 
masses and distribution amplitudes. The definitions of the hard function h e and the hard scales t\' 2 can also be found 
in Ref. [62]. Note that the factors S t in the second terms depend only on x B , instead of x B — 77 as in Eq. (1A11D . 
because we considered the small 77 region for the B s —> 7r + 7r _ /i + /i _ decay here. 
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